The crossover from vibrational to rotational collectivity in heavy nuclei in the 

shell-model Monte Carlo approach 
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We use the shell-model Monte Carlo approach to study microscopically the crossover from vi- 
brational to rotational collectivity in families of even-even samarium and neodymium isotopes. In 
particular, we identify a signature of this crossover in the low-temperature behavior of (J 2 )t, where 
J is the total spin and T is the temperature, and find it in agreement with its values inferred from 
experimental data. We calculate the state densities and find them to agree very well with experimen- 
tal data. We also calculate a collective enhancement factor from the ratio of the total state density 
to the intrinsic state density as calculated in the finite-temperature Hartree-Fock-Bogoliubov ap- 
proximation. The decay of this enhancement factor with excitation energy is shown to be correlated 
with the pairing and shape phase transitions in these nuclei. 

PACS numbers: 21.60.Cs, 21.10.Ma, 21.60.Ka, 27.60. +j, 27.70.+q 



Introduction. The shell-model Monte Carlo (SMMC) 
method 0, Q has been used successfully for the micro- 
scopic calculation of statistical and collective properties 
of atomic nuclei, and in particular level densities 
The SMMC approach enables for fully correlated shell 
model calculations in much larger configuration spaces 
than those that can be treated by conventional diag- 
onalization methods. Using a proton-neutron formal- 
ism and stabilization techniques in the canonical ensem- 
ble, SMMC has been applied successfully to the well- 
deformed rare-earth nucleus 162 Dy Here we apply 
the SMMC approach to study microscopically families of 
even-even rare-earth isotopes. Such isotopic families are 
known to exhibit a crossover from vibrational to rota- 
tional collectivity as the number of neutrons increases 
from shell closure towards the mid-shell region. This 
crossover corresponds, in the thermodynamic limit, to 
a quantum phase transition from spherical to deformed 
nuclei. However, the microscopic description of such a 
crossover in the framework of a truncated spherical shell 
model has remained a major challenge. The dimension- 
ality of the many-particle shell model space required to 
describe heavy rare-earth nuclei is many orders of magni- 
tude beyond the capability of conventional diagonaliza- 
tion methods. The SMMC approach, while capable of 
treating such large model spaces, does not provide the 
detailed spectroscopy that is used to identify the appro- 
priate type of nuclear collectivity. 

Here we show that the crossover from vibrational to 
rotational collectivity can be identified through the tem- 
perature dependence of (J 2 )t, with J being the total 
nuclear spin and T the temperature. This thermal ob- 
servable can be calculated in the SMMC method and we 
use it to demonstrate that the above crossover can be de- 
scribed microscopically in the framework of a truncated 
spherical shell model approach. Furthermore, we find 



that the temperature dependence of (J 2 )t agrees well 
with its values extracted from experimental data. We 
also calculate the total state densities for the correspond- 
ing samarium and neodymium isotopes and find them in 
very good agreement with experimental state densities. 

Vibrational and rotational collective states account for 
a significant fraction of the total state density up to mod- 
erate excitation energies. Both empirical models and 
combinatorial models of level densities often utilize phe- 
nomenological enhancement factors to account for such 
vibrational and rotational collectivity 1(1 11|- Although 
various expressions for collective enhancement factors are 
available in the literature [12, El , it is highly desirable 
to study such enhancement factors microscopically. In 
particular, little is known about the decay of collectivity 
with excitation energy although it plays an important 
role in fission reactions [l3|. Here we define a total col- 
lective enhancement factor as the ratio between the total 
state density to the intrinsic state density as obtained 
within the thermal Hartree-Fock-Bogoliubov (HFB) ap- 
proximation and study microscopically its decay with ex- 
citation energy. We find that the damping of the vibra- 
tional and rotational collectivity with excitation energy is 
correlated, respectively, with the pairing and shape phase 
transitions in these nuclei. 

Model space and interaction. Here we use the same 
single-particle model space as in Ref. Q, namely O57/2, 
W5/2, 14/2, 2si/ 2 , O/in/a, I/7/2 for protons, and 0h n /2, 
O/19/2, I/7/2, I/5/2, 2p 3/2 , 2p 1/2 , 0i 13/2 , and lg 9/2 for 
neutrons. This model space is larger than one major 
shell for both protons and neutrons, and was determined 
by examining the occupation probabilities of spherical 
orbitals for well-deformed rare-earth nuclei [9j . 

The single-particle energies in the configuration- 
interaction shell-model Hamiltonian are determined so 
as to reproduce the single-particle energies of a spherical 



Woods-Saxon plus spin-orbit potential in the spherical 
Hartree-Fock approximation. The effective interaction 
consists of monopole pairing and multipole-multipolc 
terms (quadrupole, octupole and hexadecupole) Q 
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normal ordering, 

m .„ are monopole pair oper- 



with dj m = (— y +m a,j— m is the mass 2 A -pole operator. 
The pairing coupling strengths are parametrized by 
9u =r T' 9u with g p = 10.9/ Z and g n = 10.9/JV (Z and iV 
are the number of protons and neutrons, respectively). 
The latter are determined so that the pairing gaps 
calculated in the number-projected BCS approxima- 
tion could reproduce the experimental even-odd mass 
differences for spherical nuclei in the mass region 
The factor 7 is an effective suppression factor of the 
overall pairing strength, part of which may be ascribed 
to the fluctuations induced by pairing correlations 
beyond the number-projected BCS approximation. The 
multipole-multipole interaction terms we include in 
([1]) are the quadrupole, octupole and hexadecupole 
terms (i.e., A = 2,3,4). Their strengths are given by 
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while 



Xx = X ' kx, where x is determined self-consistently [14 1 
and kx are renormalization factors accounting for core 
polarization effects. 

In general, the moment of inertia I of the ground-state 
band for a deformed nucleus is sensitive to 7 .while the 
slope of lnp(E x ) is sensitive to k 2 [9(. In Ref. [9j we have 
adopted the values 7 — 0.77, k 2 — 2.12, k 3 — 1.5 and 
k 4 = 1 for 162 Dy. We have studied families of samarium 
( 148_155 Sm) and neodymium ( 143 ~ 152 Nd) isotopes (both 
even and odd) and found that a more appropriate choice 
to reproduce the overall experimental systematics is kz = 
1, while 7 and k 2 are parametrized by a weak and smooth 
^-dependence 7 = 0.72 - 0.5/[(iV - 90) 2 + 5.3] and k 2 = 
2.15 + 0.0025(7V-87) 2 . 

The crossover from vibrational to rotational collectiv- 
ity. At low temperatures, the observable (J 2 )t is domi- 
nated by the ground-state band. Assuming a vibrational 
or rotational ground-state band with an excitation en- 
ergy E 2 + of the first excited J = 2 + , we find [^, [l5| 



(J 2 )t 



30 



7 — - — 1 + 2 vibrational band 



E. 



2+ 



-T 



rotational band 



Thus, the low-temperature behavior of (J 2 )t is sensitive 
to the type of collectivity and can be used to distinguish 
between vibrational and rotational nuclei. 

In Fig. [TJ we show the low-temperature behavior of 
(J 2 )t for the even-even samarium isotopes 148_154 Sm. 
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FIG. 1: (J 2 )t as a function of temperature in a family of 
even-cven samarium isotopes 148 - 154 Sm. The SMMC results 
(open circles) are compared with the experimental results de- 
duced from known low-lying levels (dashed lines) and from 
the additional contribution of higher levels described by an 
experimental BBF level density (solid lines). 



While 148 Sm exhibit a rather soft response in (J 2 )t, typi- 
cal of vibrational nuclei, this response becomes gradually 
more enhanced and rigid as we proceed to the heavier 
samarium isotopes. In 154 Sm, (J 2 )t acquires the char- 
acteristic linear behavior of a rotational nucleus . In- 
deed, the vibrational band formula in ([2]) can be well 
fitted to the SMMC results for (J 2 ) T in 148 Sm with 
E^ = 0.538 ± 0.031 MeV, in good agreement with the 
experimental value of E^. p = 0.550 MeV. Fitting the ro- 
tational band formula in © to the SMMC results for 
154 Sm, wc find = .087 ± 0.006 MeV, in agreement 
with the experimental value of E^_ p = 0.082 MeV. 

At sufficiently low temperatures, we can extract the 
experimental values of (J 2 )t from 
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where the summations are carried over the experimen- 
tally known low-lying energy levels i with excitation en- 
ergy Ei and spin Ji . These are shown by the dashed lines 
in Fig. [T] However, since the experimental level scheme 
is incomplete above a certain energy, Eq. ^ underesti- 
mates the correct experimental value of (J 2 )t above a 
certain temperature. We can obtain a more realistic es- 
timate by using the discrete sum over energy levels up 
to a certain energy threshold En (below which the ex- 
perimental spectrum is complete), and estimate the con- 
tribution of levels above En in terms of an average level 



density p(E x ) that is parametrized with the help of avail- 
able experimental data. We then have 

1 / N 

(J 2 )t = ^^J 4 (J l + l)(2J l + l) e - £ */ T + 

H dE xP (E x )(J 2 ) E „e- E */ T ) , (4) 

J E N J 

with Z{T) = Ef (2J i + l)e- E */ T + f~ dE x p(E x )e~ E */ T 
is the corresponding experimental partition function. 
Here (J 2 )_e x is the average value of J 2 at a given excita- 
tion energy E x . For the level density we use a backshifted 
Bethe formula (BBF) with single-particle level density 
parameter a and backshift parameter A, extracted from 
the neutron resonance data (when available) and count- 
ing data at low energies. Using the spin-cutoff model 
(obtained assuming random coupling of the individual 
nucleon spins [16|), we have (J 2 )_e x = 3(J 2 )_e x = 3a 2 (E x ) 
where a 2 is the spin-cutoff parameter. The latter is esti- 
mated from a 2 = IT/fr 2 using T = [(E — A) /a} 1 / 2 and 
a rigid-body moment of inertia 1 w 0.015 A 5 ^ 3 h 2 . The 
corresponding results for (J 2 )t are shown by the solid 
lines in Fig. [T] We find reasonable good agreement of 
these experimental estimates of (J 2 )t with the SMMC 
results along the crossover from 148 Sm to 154 Sm. This is 
a confirmation that our truncated spherical shell model 
Hamiltonian can describe the crossover from vibrational 
collectivity in 148 Sm to rotational collectivity in 154 Sm. 

The determination of the ground-state energy for even- 
even isotopes. An accurate estimate of the ground-state 
energy Eq is crucial in obtaining the excitation energy 
E x = E — Eq necessary for the calculations of level den- 
sities. Because of the low excitation energies in the heavy 
rare-earth nuclei, we have carried out calculation of the 
thermal energy up to an inverse temperature value of 
(3(= l/T) ~ 20 MeV -1 9]. The ground-state energy 
can then be reached by extrapolating the thermal energy 
to the limit j3 — oo. In the vibrational and rotational 
limits we have used expressions for the low-temperature 
energy in the ground-state band approximation 0, [l5| 
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We have extracted the ground-state energy Eq of 
148 Sm and 154 Sm from a one-parameter fit to the low- 
temperature SMMC values of E(T) using, respectively, 
the vibrational and rotational band formulas in ([5]) [in 
the vibrational case, E 2 + was already determined from 
the fit of (J 2 )t using Eq. (J2J)]. For other nuclei in the 
crossover between the vibrational and rotational limits, 
we have estimated Eq by taking average value of E(T) 
at sufficiently low temperatures. We have found Eq = 
-235.690±0.015, -255.772±0.019, -275.849±0.018 and 
-295.398 ± 0.024 MeV for 148 Sm, 150 Sm, 152 Sm, and 
154 Sm respectively. 
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FIG. 2: Total state densities in the even-even 148_154 Sm iso- 
topes. The SMMC results (open circles) are compared with 
level counting data (histograms), neutron resonance data (tri- 
angles) and the BBF parametrization of the experimental 
data (solid lines). Also shown are the HFB level densities 
(dashed lines). The neutron and proton pairing transitions 
are indicated by arrows and the shape transition by thick ar- 



State densities in samarium and neodymium isotopes: 
theory and experiment. In Fig.[2]we show the total state 
densities as a function of the excitation energy E x for the 
even-even samarium isotopes ( 148_154 Sm). The SMMC 
state densities (circles), calculated using the method of 
Ref. are compared with experimental data that con- 
sist of level counting data at low excitation energies (his- 
tograms) and, when available, the neutron resonance 
data at the neutron threshold energy (triangles). For nu- 
clei with neutron resonance data, we have also included 
a BBF state density [13] (solid lines) whose parameters 
a and A are determined from the level counting and the 
neutron resonance data [Hj]. 

In the SMMC calculations, the imaginary time (3 is 
discretized in steps A/3. To take into account systematic 
errors due to discretization, we have done calculations 
for two values of A/3 (1/32 and 1/64 MeV^ 1 ). For (3 < 3 
MeV -1 the energies are extrapolated linearly to A/3 = 0, 
while for larger values of (3 an average value is taken Q . 

For comparison, we also show in Fig.[2]the level density 
calculated from the finite-temperature HFB approxima- 
tion (dashed lines) with the same Hamiltonian. The HFB 
level density accounts only for intrinsic states, and there- 
fore the enhancement observed in the SMMC state den- 
sity originates in vibrational and rotational bands that 
are built on top of these intrinsic states. The kinks in 
the HFB level density are associated with the proton and 
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FIG. 4: Total collective enhancement factor K (see text) in 
the even-even 148-154 Sm isotopes as a function of excitation 
energy E x . Arrows are as in Fig. [21 
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FIG. 3: Total state densities in the even-even 144 - 152 Nd iso- 
topes. Results shown are as in Fig. [2] 

neutron pairing phase transitions (arrows) and the shape 
phase transition (thick arrows). 148 Sm is spherical in 
its ground state and undergoes pairing transitions only. 
150 Sm has a non- vanishing deformation in its ground 
state and undergoes also a shape transition to a spherical 
shape at E x re 12.5 MeV. The ground-state deformation 
continues to increase with mass number in 152 Sm and 
154 Sm, and the shape transitions occur at higher excita- 
tion energies (outside the energy interval shown in the 
figure) . 

We have carried similar studies of state densities for 
the family of even-even neodymium isotopes ( 144 ~ 152 Nd). 
The results, shown in Fig. [31 exhibit similar trends to 
what we have observed in the samarium isotopes. 

Collective enhancement factors. We define a collective 
enhancement factor by the ratio K = psmmc/phfb- In 
Fig. [4] we show K (on a logarithmic scale) versus excita- 
tion energy E x for the same samarium isotopes of Fig. [3J 
but up to higher excitation energies of E x ~ 30 MeV to 
view effects of the shape transition. 

148 Sm is spherical in its ground state and the observed 
collective enhancement must be due to vibrational col- 
lectivity. This collectivity disappears (i.e, K re 1) above 
the proton pairing transition. The other samarium iso- 
topes shown in Fig. 2] are deformed in their ground state 



and K exhibits a local minimum above the pairing tran- 
sitions, which we interpret as the decay of the vibrational 
collectivity. The rapid increase of K above the pairing 
transitions originates in rotational collectivity. This col- 
lectivity reaches a plateau as a function of excitation en- 
ergy and then decay gradually to K ~ 1 in the vicinity of 
the shape transition (thick arrow) when the nucleus be- 
comes spherical and no longer supports rotational bands. 

Conclusions. We have carried out SMMC calculations 
for isotopic families of the even-even rare-earth nuclei 

148-154 Sm and 144-152 Nd Ugmg the observab l e (j2} Tj 

whose low-temperature behavior is sensitive to the spe- 
cific type of nuclear collectivity, we have demonstrated 
that a truncated spherical shell model can describe the 
crossover from vibrational to rotational collectivity in 
heavy nuclei. We have also calculated the total SMMC 
state densities and found them to be in very good agree- 
ment with experimental data. We have extracted micro- 
scopically a collective enhancement factor defined as the 
ratio between the SMMC and HFB state densities. The 
damping of vibrational and rotational collectivity seems 
to be associated with the pairing and shape phase tran- 
sitions, respectively. 
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